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^ Abstract. Wc study the heat content asymptotics with either Dirichlet or 

Robin boundary conditions where the initial temperature exhibits radial blowup 
near the boundary. We show that there is a complete small-time asymptotic 
, expansion and give explicit geometrical formulas for the first few terms in the 

expansion. 



in 



1. Introduction 



1.1. The heat content. Let M be a compact Riemannian manifold of dimension 
m with smooth boundary dM, let dx and dy be the Riemannian measures on M and 
■ on dM, respectively, and let D be an operator of Laplace type on a smooth vector 

bundle V over M. To impose suitable boundary conditions, we define the Dirichlet 
boundary operator Bx>4> '■= 4>\dM- The operator D defines a natural connection V 
as we shall discuss presently in Lemma 11.11 Let S be an auxiliary endomorphism 
of the vector bundle V\qm- Let 

>- ■ B-]i4> '■= {4>;m + S<j))\dM 

be the Robin boundary operator where (f> :m denotes the covariant derivative of <j) 
\q \ with respect to the inward unit normal vector field. Let B be either the Dirichlet 

C^) • or the Robin boundary operator; the associated boundary conditions are defined by 

\ setting Bcj) = 0. It is well known that the heat equation 

(fit + D)u B (x; t) =0, Bu(-,t) = 0, \imu B {-; t) = </>(■), 
00 . '1° 

has a unique classical solution for a wide class of initial temperature distributions 

</>. We set u — e~ tDB <fi where Db is the associated realization of D. The operator 

e~ tDs has a kernel ps{x,x;t) which is smooth in (x,x;t) such that 

C$ ■ u B (x;t) = I p B (x,x;t)4>(x)dx . 

If, for example, D = A := 5d is the scalar Laplacian, then one may take a complete 
spectral resolution {A;,0i} of A B and express 



Pb 



(x,x;t) = tXi (j)i{x)(j>i{x) . 



Let (•,•) denote the natural pairing between V and the dual bundle V. The 
specific heat p of the manifold is a smooth section of V and the heat content /3 is 
given by: 

f}((/),p,D,B)(t) := / (u B (x;t), p(x))dx = / / (p B {x, x; t)(j>{x), p(x))dxdx . 
Jm Jm Jm 

Although in most practical applications it is customary to take V and V to be the 
trivial line bundle and D — A, it is necessary to work in this greater generality as 
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we shall see presently in Section [3] If dM is empty, we shall omit the boundary 
condition B from the notation as it plays no role. To simplify the notation, we shall 
write Db, Pb, /?(-, -,D,B), and ub for the most part except where it is useful to 
emphasize which boundary condition appears. 

1.2. Geometric preliminaries. The following formalism will enable us to work 
in a tensorial and coordinate free fashion. We adopt the Einstein convention and 
sum over repeated indices. Choose a system of local coordinates x = (xx, ...,x m ) 
for M and choose a local trivialization of V. Let g M „ := g(d Xft , d x J) and let g^ v be 
the inverse matrix. As D is of Laplace type, there are matrices A" and A® so that: 

(l.a) D = -{g' iV lAd Xfl d Xv +A v 1 d Xu +A } . 

If V is a connection on V, we use V and the Levi-Civita connection to covariantly 
differentiate tensors of all types and let ';' denote multiple covariant differentiation. 
We let 4>. tllv be the components of V 2 0. If E is an auxiliary endomorphism of V, 
we define the associated modified Bochner Laplacian: 

D{g,V,E)<l>:=-gr<l>., vll -E(l>. 

Let r and r M „ CT be the Christoffel symbols. We then have [9]: 

Lemma 1.1. If D is an operator of Laplace type, then there exists a unique con- 
nection V on V and a unique endomorphism E on V so that D = D(g, V, E). The 
connection I- form uofV and the endomorphism E are given by: 

(1) ^ = i( 9 ^ + ff -r^id). 

(2) E = A Q - g^(d Xv ou^ + uj^uj v - lo.T^). 

We use the dual connection to covariantly differentiate the specific heat p; note 
that the connection 1 form H> u for V is the dual of —uj v . Thus 

(l.b) Vau„ =d Xtl -i(g^ 1 +g^T (7Sfl id) and Dp = - (g^p. tfiu + Ep) . 

Near the boundary, choose an orthonormal frame {ei,...,e m } for the tangent 
bundle so that e m is the inward unit geodesic normal; let indices a,b range from 1 
to m — 1 and index the induced orthonormal frame {ex, ...,e m _i} for the tangent 
bundle of the boundary. We let ':' denote the components of tangential covariant 
differentiation defined by V and the Levi-Civita connection of the boundary. Let 
Lab '■= <?(V 'e a e&,e TO ) = r a (, m be the components of the second fundamental form. 
The difference between ';' and ':' is then measured by L. For example, the following 
relation will prove useful subsequently: 

Q \ D4> = - {4>:aa + <fi;rmn ~ L aa (j). m + E(j>), 

Dp = —{p:aa+ P;mm — L aa P;m+ Ep) . 

Let Ric denote the Ricci tensor. Let B^ be the dual Robin boundary operator; it 
is defined by the dual connection V and dual endomorphism S. 

1.3. Heat content asymptotics in the smooth setting. One can use the cal- 
culus of pseudo-differential operators developed in [Till HH 03] *° snow that: 

Theorem 1.2. Let <\> 6 C°°(V) and let p € C 00 ^). There is a complete asymptotic 
expansion as t I of the form: 

00 (—f\ n ( °° ( 

(3(^p,D,B)(t)~Y,—r (^D n p)dx + Y / t (1+k)/2 [3 9 k M {^ P ,D,B)dy. 

n=0 U - Jm k=0 JdM 

The heat content coefficients /3® M are locally computable and arc given by local 
geometric invariants. 

Theorem 1.3. Let <j) E C°°(V) and let p e C°°(V). 

(1) With Dirichlet boundary conditions, one has that: 
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( a ) IdM /3o M (^ P, D, B v )dy = f gu {<l>, p)dy. 

( b ) J dM P! M (ct>,P,D,Bv)dy = -J aM {(<t>;m,p) ~ \{L aa <j>,p)}dy. 

( C ) I aM P2 M {<t>,P,D,Bv)dy = -^/ aM {(f<^mm,p) + f<0,P ; mm> 
-|£oa(0, P);m + P) ~ (0:a, P:a) + ((^L aa L bb 

-\LabLab - \ Ric mm )0, p)}dy. 
(2) H^i/i Robin boundary conditions, one has that: 

(a) J dM $ M (4,,p,D,B n )dy = 0. 

( b ) Jom Pi M (<t>, P» £> » = Jg M <0, B K p)dy. 

(c) J aM & 9M (0, P, D, B n )dy = ^ J 9M (B R ^ B n p)dy. 

With Dirichlct boundary conditions, /3f M is known and partial information con- 
cerning [3f M is available. With Neumann boundary conditions, fa, fa, and fa are 
known. We refer to [9] for further details; our indexing convention here is slightly 
different from that employed in [9]. 

1.4. Singular initial temperatures. We refer to [2] for earlier work in the sin- 
gular setting as the results of that paper provide the motivation and the starting 
point for this present work. We begin by using the geodesic flow defined by the unit 
inward normal vector field to define a diffeomorphism for some e > between the 
collar C e :— dM x [0, e] and a neighborhood of the boundary in M which identifies 
dM x {0} with dM ; the curves r — > (yo,r) for r 6 [0, e] are then unit speed geodesies 
perpendicular to the boundary and r is the geodesic distance to the boundary. We 
fix a smooth cutoff function x — x( r ) on Ce so that x = 1 near r — and so that 
X = near r = e. 

Let V be a connection on a bundle W over C e . Let G C°°{W\qm)- We use 
parallel translation along the normal geodesic rays to extend ip to a section of W 
over C e . We shall denote this extension by tp(y) to emphasize the fact that Vg r tp = 
on C e . If W = V, we use the connection V defined by D; if W = V, we use the 
dual connection V defined by D. We refer to [9] for further details. 

Fix a e C. Let be a smooth section to V on the interior of M such that 
4>r a S C°°(C e ); the parameter a controls the growth (if Re(a) > 0) or decay (if 
Re(a) < 0) of <f> near the boundary assuming that <fir a does not vanish identically 
on the boundary. We assume the specific heat is smooth so p 6 C°°(V). We may 
then expand <p and p on C e in the form: 

oo oo 

(l.d) 4>(y,r) ~ 5Z<My)r l ~ a and p(y, r) - as rlO. 

i=0 i=0 

The coefficients cf>i and pi are then uniquely specified by the requirement that 
Vd r (/>i — and Vgi r p, = 0. If D = A, then the associated connection is flat and, if 
a = 0, the expansion of Equation (|l.d[) is just the usual Taylor series expansion of 
the functions <p and p. In particular, 

00 = (r Q 0)U/, 01 = {V ar (r Q 0)}| aM , 02 - H(^9r) 2 (r a m\dM 

P0 = p\dM, Pl={^drp}\dM, P2 = \{{^drYp}\dM ■ 

Fix t > 0. Let a; e M and let x = (y, r) £ C e . Suppose first that B = Bx> defines 
Dirichlet boundary conditions. Then Pb{x, (y, r), i)|f=o = 0. Since ps is smooth 
for t > and C e is compact, we may use the Taylor series expansion of ps to derive 
the estimate: 

\p B {x,{f,y);t)\<C{t)f on C e . 
If Re(a) < 2, then the integral us{x]t) — f M Pb(x, x;t)(j)(x)dx is convergent and 
bounded in x. Consequently the heat content 

fa<j>,p,D,B){t) := f (u B {x;t),p(x))dx 
Jm 
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is well defined for t > 0. If Re(a) < 1, then G L 1 and the initial heat content 
(3(4>, p, D, B)(0) = f M ((j),p)dx is finite. If, however, 1 < Re(a), then this integral 
may be divergent and the initial heat content can be infinite. Still, sufficient cooling 
near the boundary takes place for ub to be in L 1 for any t > 0. A similar phe- 
nomenon occurs in the setting of non-compact Riemannian manifolds with infinite 
volume and with regular boundary and initial temperature <f> = 1 [HE]. 

As this cooling phenomenon does not occur with Robin boundary conditions 
B = B-k, we shall always assume Re(a) < 1 in this instance. 

It is important to observe that although we are primarily interested in positive 
real a, it is necessary to consider complex values of a to justify some analytic 
continuation arguments. It is also necessary to permit Re(a) < to justify some 
computations in Sections [3] and 01 these values are of interest in their own right 
since <fi is not smooth if a is not an integer. 

If 1 < Re(a) < 2, we must regularize the integral J M (<p, p)dx since the integral 
may be divergent. The Riemannian measure is not in general product near the 
boundary. Since, however, dx = dydr on the boundary of M, we may decompose 

(0, p)dx = (0 O , p )r- a dydr + 0(r l - a ) . 

For Re(a) < 2, define: 

(l.e) 2Re g (0, p) ■= (4>,p)dx+ {(</>, p)dx - (0o, po)r~ a dydr} 

Jm-c c Jc, 



(<t> ,po)dy x 

dM 



% if a + 1, 

ln(e) if a = 1 . 



This is clearly independent of e and agrees with f M (<fi,p) if Re(a) < 1. Briefly, the 
regularization 2R Cg (0, p) is a meromorphic function of alpha with a simple pole at 
a = 1. When a = 1, then 2R Cg (</>, p) is defined as the constant term in the Laurent 
expansion at a = 1, thus dropping the pole. 

The following is the main analytic result of this paper: 

Theorem 1.4. If B = Bx>, assume Re(a) < 2; if B = Br, assume Re{a) < 1. 
(1) If a =^ 1, then there exists a full asymptotic expansion as t J, of the form: 

0(<f>,p,D,B)(t) ~ f^t^l lRcg ^D n p) 

71 = 



OO „ 

k=o JdM 



(2) If a = 1 (and hence B = Bx>), then there exists a full asymptotic expansion 
as 1 10 of the form: 

/3(0, p, D, B)(t) ~ V { —^I Rcg (^ D n p) 

n=0 

E ^ f P> D > B ) + ln W/3f M (^ P, D, B)}dy. 

t_n JdM 



k=0 



(3) There are natural tangential bilinear differential operators fi^aij' which 
are holomorphic for a/1, so that 

[3™ (</>, p, D,B)= E ffiXui Pi ,D,B). 



i+j<k 



(4) 



I DM 



f3 d k M (cf>, p, D, B)dy = { Um^o, (D n p) )dy if k = 2n 

I if k is odd. 
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We note that B-jip = pi+Sp$. Theorem l 1 .31 generalizes to this setting to become: 
Theorem 1.5. Set c a := 2 1 " < T * 

(1) //a^l, i/Re(a) < 2, and if B = B v , then: 

( a ) JaM $>,a f (^ P> A = c Q f dM {<l>o, Po)dy. 

( b ) Jam Pi^^i P> D > B v)dy = c a _ x $ dM {<t>i - \L aa <j) , p )dy. 

~ 2(a- Q l7( 3 a-2) ( E< t>0,PQ) + ( Q _i) 2 ( Q _ 2 )(^0,P2> - ( a _l)(a,_2) { L aa<t>0, Pi) 

+ 2(a-°l)"(a-2) (^0:oi P0:a) - 4(a-2) ( Ric ™™ 00, A)) 

+ 8(a"-2) (LaaLbb&O, Po) ~ 4(^-2) (LabL a b4>0, Po)}dy. 

(2) Let 7 6e Euler's constant. If a — 1 cmd if B = Bx>, then 

0(4>, P, D,B v )(t)~ 2R, cg 0, p) - | ln(t) f gM {<f>o,po)dy + f dM %{(fo, Po)dy 
+ il/2 /ai;{-^(^./ , o) + ^{Laa<t>Q,Po)}dy 
+t{-J Reg (0 ! Dp) + i ln(i) J dM {^o, (Dp) )dy} 

+*/9M{i^ ' _ (^)°) _ (^2,Po) + \ {L aa <t>l, Pq) + (00, P2> 

-§<Wo,Pi»dy +o(i 3 / 2 ). 

(3) Ifa^Q, if(Re)(a) <1, and if B = B n , then: 

( a ) J 9M P™(<t>,P,D,B n )dy = 0. 

( b ) JaM /ta P> A - ^c Q+1 J aM (0 o , B KP )dy. 

(c) J dM (3™(<f>,p,D,B K )dy 

= 3^ c " Jsm^ 1 ~ "O^ 1 + - %L aa (j> ,B n p)dy. 

Remark 1.6. We note that setting a = in Assertion (1) of Theorem 11.51 yields 
Assertion (1) of Theorem 11.31 and that taking the limit as a — > in Assertion (3) 
of Theorem 11.51 yields Assertion (2) of Theorem 1 1.31 

To avoid subscripts on subscripts, we shall for the most part simply refer to Ms, 
Db, and ps when no danger of confusion is likely to ensue; however, we shall use 
the notation Db^ and Db k in Section [4] when we must deal with two different 
boundary conditions. 

I. 5. Outline of the paper. In Section[2l we use the calculus of pseudo-differential 
operators to establish Theorem 11.41 (1-3); we postpone the proof of (4) as it will 
follow as a scholium to the proof of Lemma l3~Tl it can also be deduced directly from 
Lemma 12.61 We shall restrict to Dirichlet boundary conditions as the analysis is 
similar for Robin boundary conditions. In Section[31 we apply invariance theory and 
the functorial method to prove Assertions (1) and (2) of Theorem II .51 an essential 
input is the calculation of a single coefficient performed in [5] and a special case 
computation on the half- line. In Section dj we establish Assertion (3) of Theorem 

II. 51 We plan in a subsequent paper [5] to undertake a similar analysis of the heat 
trace asymptotics with a singular smearing function. 

The special case a = 1 in Theorem 11.51 can be derived from the case a =/= 1 
by a straightforward but tedious analytic continuation; the poles at a = 1 which 
arise from the terms in Assertion 1 of Theorem 11.51 are cancelled by those involved 
in the regularization lR eg (0, D n p). The derivation of the t° and t°\n(t) terms in 
this way is sketched in Section [2~T1 the other terms use Equations (|l.b|) . (|l.c[) . and 
(|3.bj) below. We will give an independent derivation of the case a = 1 to show how 
scaling arguments can be applied when logarithmic terms are present. 
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2. Pseudo-Differential Operators 

We suppose throughout Section [5] that Re(cc) < 2 and that B = Bx>. We use the 
identity 



{e- tUB (j),p)dx = / {<j),e- tUB p)dx 

I M JM 

to interchange the roles of tj> and p. This causes no difficulty as 

p(cl),p,D,B v )(t)=/3{p,ci>,D,Bv)(t). 

We shall assume throughout this section that <j> is smooth and that 

P(y,r) = Po(y)x( r ) r ~ a on C e . 

The more general case where p ~ (po(y) + Pi(y) r + •■■)x( r ) r Q then follows simi- 
larly; if <p and p vanish to high order on the boundary, the corresponding boundary 
contributions vanish to high order in t. Frequently in this section, we will let (3 be 
a multi- index rather than denoting the total heat content; we shall also let 7 both 
be a multi-index and Euler's constant. We apologize in advance for any confusion 
this may cause. We let lR Cg be the regularization defined in Equation (| 1 .e|) where 
we interchange the roles of <j> and p. 

Our fundamental analytical result is the following: 

Theorem 2.1. Adopt the notation established above. There are differential op- 
erators Gfc of order at most k defined on M and tangential differential operators 
Bfe(a) and of order at most k defined on dM such that in closed subsectors of 
{p € C : I arg(/i)| < ^tt} one has: 

(1) Let a ^ 1 and let Re(a) < 2. Then Bo (a) = — T(l — a), Bfc(a) is holomor- 
phic in a, and 

00 00 „ 

{(Db+p 2 )- 1 ^^) ~^~ 2 ~ fc WG fe <M + 5> Q ~ 3 ~ fe / (®k(a)<j>,p )dy. 
fc=0 k=0 JdM 

(2) Ifa = l, then: 

((Db + pY'^p) - Yu- 2 - k {iR Cg (G k 4>,p)+ I (n k (l)4>, Po )dy 



/ J I 

k=0 



OM 



Tn(/z) / (h k 4>,p )dy 

JdM 



Here is a brief outline of the proof. In Section [2711 we study the half- line M.+ . In 
Section [2~2l we study the half-space K™. In Section [2~3l we complete the proof by 
considering the case of manifolds. Theorem II .41 then follows. The interior integrals 
are evaluated using previous techniques. Furthermore, for a^l, one has 



(0, p, D, B V ) = r(^f±5) (B k (a)h p ) . 

In particular, we may use the duplication formula for the Gamma function to es- 
tablish Assertion (la) in Theorem 1 1.51 by computing: 



(3™(<j>,p,D,Bv) = -§^(h,Po) = 2 x -«r (2=s) 7 =A- Ty (0 o , Po ) . 

2.1. The half line. The case of gives the basic outline of the proof in the 
general case. Let 77 be the Fourier transform variable related to r. The Fourier 
inversion formula then becomes: 

POO />oo 

,£(77) = / e~' /=Trri 4>{r)dr and (j>(r) = / e^ 1 ^ 4>(<r))dr) 
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where we set dr/ :— drj/2'K. For D — —d 2 , we may then write the Dirichlet resolvent 
as a pseudo-differential operator in the form: 



{D B +^r^{r) 



e v lr ^ C - 2 {ri,^)<i'{'n)d'q+ I d- 2 (r,r),[j,)(f>(ri)dr) 

J — CO 

= : [Op(c_ 2 )+Op'(d_ 2 )Mr). 

Here c_2 = (^ 2 +A* 2 ) _1 j and cL 2 (r, r], p) = — (rj 2 +p 2 )~ 1 e~ fJ,r is the bounded solution 
of 

(-d 2 + p 2 )d- 2 = and d_ 2 | r =o — —C-%{r}, p) ■ 
The kernel k(r, s, p) of (As + p 2 )^ 1 is thus, for r > and s > 0, given by 

1 r 



k(r,s,p) 



2fi 



e \r-s\fJ. _ e ~(r+s)n 



We now consider the heat conduction problem: 

<fi(r) = 1 and p(r) = x( r ) r ~ Q i where x is smooth and 

1, < r < 1/3, 
0, 2/3 < r. 

Suppose that Re (a) < 1. One then has that: 



(2.a) 



X(r) 
. One 

// J ((£> B +M 2 )"V.P)=i" 2 / \"">'- ■" / l<-(r.,.,,)d,dr 

X {r)r- a (l-e- rll )dr 

^ a dr + J (x(r) - l)r~ a dr - / / "',lr 
( X (r) - l)e- rll r- a dr 



+ / (X(r) - l)r- a dr 



/i a - x r(l - a) + 0(ji-°°) . 



XI- a) 

The term in bracket gives lR g (0, p) for all a ^ 1 , hence 

(2.b) p 2 ((D B + M 2 )"V, p) = 2Reg(0, P) - M a_1 r(l - a) + O^" 00 ), a ^ 1 

As a — > 1 we note that T(l — a) = — 7 + ... where 7 is Euler's constant. 
Consequently 

1/(1 - a) - /i a - x r(l - a) -» ln(/x) + 7 . 
So as a — > 1, Equation (|2.ap gives 

(2.c) M 2 ((D B + M 2 )-V,P) =ln(M)+7+^Rc g (0,p) + O( A1 -°°), a = 1. 
We pass to the heat content by a contour integral 

1 



27TV^T 



e tx ({D B + X)- 1 cj) 1 p)dX, 



where C is a "keyhole contour" consisting of two rays {re 



±V=T(7T-e) 



r > R} and a 



circular arc {Re 



\6\ < 7r - e}. We then get 



= / - rg^) + Q(*°°), (a^l) 



(2.d) (3{cj),p,-d 2 ,B D )(t) = , , - , 

[l^y-JhW + ^ + Off), (a = l) 

This last formula is valid also for the same functions </> and p on the interval [0, 1] 
with Dirichlet conditions at both ends; since p = near r = 1, the boundary 
correction from r = 1 is 0(t°°). The constant appearing in (|2.dj) for a / 1 jibes 
with the constant c Q defined in Theorem 11.51 by Legendre's duplication formula. 
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Thus, for this special case, we have established Assertions (1) and (2) of Theorem 
PI 

2.2. The half-space M™. We use the resolvent construction described in [T3]. Let 
£ = •••! Cm-i) be the Fourier transform variables dual to y = (yi, y m -i). 
Note that g(d Vi ,d r ) = and g(d r ,d r ) = 1. We set 

9 2 (^,0 := .9(£ady a ,6d2/ b ) = S a6 (a0&6 • 
We adopt the notation of Equation (|l.a[> . The symbol of Z? is given by 

a(D) = cr 2 + u\ + cro, where <r 2 = {<7 2 (^, + ?7 2 } id, 

ax = -V^iAlta + Af 77}, and ao = -4o ■ 
The resolvent parametrix is the sum of an interior part and a boundary correction. 
Let y ■ £ = yi£,i + ... + y m — i£m— 1< The interior part of the parametrix is a finite 
sum Op(c_ 2 ) + •■• + Op(c_jv) where 

(2.e) [Op(c)4>](y,r) = J J e^^+^c^, 7-^,77)^,77)^. 

The leading term in the interior of M is 

c-2(y,r,^,r),fi) = (<j 2 + /i 2 ) -1 

and successive terms c_3,... are defined by the usual pseudo-differential calculus 
with parameter /i; see, for example, Lemma 1.7.2 [8]. In particular, cj is homoge- 
neous of degree j in the variables (£, 77, /i) and, for | arg(/j)| < jir — e, we have the 
estimate: 

(2i) I^A^il ^ const^^KI^I + \ V \ + | M |)i-lTl . 

The boundary part of the parametrix is a finite sum of operators 

which are chosen so that 

(1) L»{Op'(d_ 2 ) + ... + Op'(d_iv)} has order 1 - N, 

(2) [Op(c J -)+Op'(^)]0(tf,O) = O. 

To ensure that this second condition is satisfied, we set: 

(2-g) dj(y,0,£,r),[x) = -Cj(y,0,£,r),ti) . 

To ensure the first condition is satisfied, we begin by setting 

[-a 2 + q 2 + /J 2 ]d_ 2 = . 

Then equation Q2.g[ ) yields 

d. 2 (y, r, 77, /*) = -{q 2 + V 2 + ^)- l e- r ^+^ . 
We may then define dj inductively by an equation of the form: 

-j 

(-a 2 + q 2 + /i 2 )^^! - Y, ajkiu0y(yyd?^djd k , 

-k=2 

where k — £ + u + \ f3\ = j + 1. The coefficients ajkiupyiy) come from the Taylor 
expansion of the coefficients of D in powers of r, see [12] for details. For some 
constant C, we have: 

(2.h) dj(y,r/t,t(i,tri,tfji) = t J dj(y,r,^,r],fi), 

(2.i) r^qfa^d,- = o(iei + M + 2 + i/xD^-i^+v^^i+w) . 

Now consider the expansion of [Op(cj)<p](y, r) for <fi G 6>(M™), i.e. assume that <j) 
has an extension in the Schwarz class iS(R m ). Set 
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4>W :=ag, ir) M,r) 

and 

Cjp(y,r,Z,r),iJ,) := / = T) l/3| ^ )T)) c i (y, r, £, 77, /i) ■ 
Lemma 2.2. As — > 00 m closed subsectors o/{/ieC: | arg(/x)| < 5^}, 
[Op(c,)0](y,r) ~ J] +^7# (/3 W), ^ 

Cjp{y,r,n) = c i/3 (y,r,0,0,/x), and 
Ajp(y,r,(j,) = - / e~^* c Cj-£(y,r,0,C,M)dC<te. 



OO •/ —OO 



Proof. In Equation (|2.e[) . consider first the integral Let 3>(y, 77) be the Fourier 
transform in r. A Taylor expansion of $(y, ?]) in powers of y ~ y gives 

/OO pOC _ 

-OO ./ — (X) 

= E ■V j7 (y ) r,0,C,M)* (7) (l/,0 

(2.j) + / / e^y-v* £ c^iy^^^^iy^mOkdy. 

\l\=K 

In this expression, 7 = (71, 7 m _i, 0) is a multi- index. The change in the order of 
integration is clearly justified if j < 1 — m. For other j, we may insert a factor of 
(1 + ICI 2 ) 1 " an d continue analytically to w = from 2 Re(u>) + j < 1 — m. 

We now multiply Equation ( |2.jp by e v/ ^" r ^ and integrate d£. The remainder 
integral is a harmless 0(|/ip +m+2 ) since Rj is bounded and for all 7' and for 
I7I = K, we have 

(y - VV' J e^y-y>^(y, r, £, m)^ = 0(|C| + lA) j ~ k+m ■ 
From the terms with I7I < K, we have 

- r e^ =rp ^ T (tf,r ) 0,C,/*)* (7) (!/ > C)^ 

OO /> CO 



1 



A Taylor expansion of 0^ in powers of s — r gives 

i /*CO pOC 

E^7/ / e^'— ^ /J (»,r,O,C,M)«*:*»0 ( ' J) (l/,r) 
plus a harmless remainder. Here /? = (71, 7 m -i, u). Writing as — J_ c 



we find the Cjp and Ajp as in Lemma 12.21 □ 
We continue our development. Let djp — {—^/—l)^d^ c dj 



Lemma 2.3. [Op' (dj)(/>](y,r) ~ ^ jrBjis^Hy^^) with 

pOC poo ____ 

Bjp{y,r,n)= I I e~^ a<: djp(y,r,0,C,(i)dCda 

JO J — oo 

Furthermore 

(2.k) [^ i/3 + S^ + C^](y,0 ) M)=0. 
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Proof. The proof is similar to the proof given in Lemma I2T21 We sketch the details 
as follows. We express 

[Op' (d j )<j>](y,r)= f f f [ e^y-^-^<d 3 (y,r^X,^(y,s)dydsd0i 



and expand 4>(y, s) in powers of (y — y, s — r). To establish Equation (|2.k() , we use 
Equation ( |2.g[ ) to see that: 

r poo 



[A jP +B jP ]{y,Q,n) = / e-^-^c 3f3 (y,0,0X^)dCd 



/.s 



oo »/ — oo 



e-^*c w (y,0 J ) C,A*K<fo 

J-oo 

= -Cjfi(y,0,0,0,fi) = -C j fi(y,0,n). 
The Lemma now follows. □ 
Next, we consider ((Db + p 2 )~ l 4>, p) as (X — > oo where p(y,r) = po(y)r~ a . Define 

/— t /-oo 
/ e-^^Cj^r^C^dCds. 
-oo J — oo 

We set := We then have 

(2.m) A jP (y,r,p) = m AJ (y, r, | M |r, //) • 

Lemma 2.4. M^e have: 

(A J/3 (y,r, A i)^«(y,r),p (y))r- Q rfr 

fc=0 ^° 
The remainder after N terms of the expansion is analytic for Re (a) < N and is 
0(\n\ 1+j - N -W) uniformly in |Re(a)| < N. 

Proof. With a change of variable |/x|r = t, we have 

'{A jP <t>W,p )r- a dr 

\a-l+j-\/3\ 



= A* 



(^(y, vImM, m')4> (/3) (v, po(2/))t- Q di . 



The Lemma follows from a Taylor expansion of Aj^(y,r,t, fjf)(f>^' (y,r) in powers 
of r. This expansion is justified as follows. In Equation (|2.1|). the integral <i£ is 
0((1 + s)~°°) and so are its derivatives in (y,r) in view of Equation (|2.fT) . Hence 
and its derivatives in (y, r) are 0((1 + □ 

For the term with Bjp<f>(P\ we define 



d~p(y,r,ti,s,ii) = / e v uc djp(y,r,£,<!;,V)d(- 

J —oo 

One then has that 

(2.n) B jf3 (y,r,fi) = dj f3 (y,r,0,s,n)ds. 

Jo 

Lemma 2.5. 

(B J/3 0(' 3 )(y,r),po(y)>r- Q dr 



^ -i n OO /-OO 

^ |/ir -i + i- fc -i^i / / {d7 p {y , rAs ^>) d ^m iy ,o) 7Po{y) y-« dsdr . 



k=Q 



JO 
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The remainder after N terms is bounded as in Lemma \2.J\ 

Proof. From Equation (|2.hp . noting that djp = (— x/^T)" 3 '^ -dj, we have 

dj (y,r/t,t£,s/t,t(i) = t 3 -^^ +1 dJ p (y,r,^s,fj,) . 
We have from Equation (|2.i|) for all 2V that 

(l + s 2 )r N dJ /3 (y,r,0,s,f/) = O(l). 
The Lemma now follows from a Taylor expansion of ^ in powers of r. □ 

Finally, since Cjp is homogeneous of degree j — \/3\ in (£, CjA 4 )) we have from 
Lemma 12.21 that: 

(C j73 (y,r, M )^(y,r),p (y)>r- a dr 



(2.o) = / (c oP {y,r,Q,Q,p')4> W {y,r),p (y))r~ a dr. 



We may add up all the expansions in Equation (|2.o|l . in Lemma l2.4i and in 
Lemma T2.5I All terms are analytic for Re(a) < 1 and we extend the expansion by 
a meromorphic continuation to Re(a) < 2 with a simple pole at a — 1. The terms 
with a singularity at a — 1 are 

poo 

lA j ~ m I (c J p(y,r,Q,0,^)cl ) ^( y ,r),p o ( y ))r- a dr 

OO 









OO pOO 



+ |^|«-i+i-|/3| / / (d^(|,,r,0,*,M / )^ ) (tf.O).A)(y)>r- a d»dr. 



'J/3 

o Jo 

We must now determine the extension to a = 1. As previously, when we regu- 
larize the interior integral, we investigate 

c i/3 (y, 0, 0, 0, //) + |/*| a_1 {^(tf, 0, 0, //) + / °° dr p (y, 0, 0, s, fi')ds} 

(2.p) ^-j ^ . 

1 — a 

From Equation (|2.mp with = 1, we have A?p(y, 0, 0, //') = Ajp(y, 0, //). From 
Equation ()2.n[) . 



dJ (y,O,O,s,p')ds = B j(3 (y,0,p r ) . 
Then from Equation (|2.kj) and Lemma l2~2l the expression in Equation (|2.p[) is 



1 — I l Q ~ 1 

^(y, o,o,o,m')— — fOT 

1 — a 

Taking the limit as a — > 1 yields 

£773(2/, 0,0,0,//) In a = l 

and consequently we have: 
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Lemma 2.6. For the parametrix P N = £^=2 [°P( c j) + °p'( d j)]j f or <t> G <S(K+), 
/or r) = po(y)r~ a , for Re(a) < 2, and /or a ^ 1, we have 

(P N(t >,p)~J2M j - l0l i I r \cMy,rM»')<P W (.y,r), Po (y))r- a drdy 
+ E ^^'^Jki [ [°° d*(AJ^\p o )(y,0 7 s,p')s k - a d S dy 

j,/3,k f ' J J ° 

+ E \^ 1+3 ^ k mu / / / ^(y,r J 0, a> /i / )r fc -°'drd a 

■(d^(y,0),p (y))dy. 

When Re(a) > 1, the divergent integrals are regularized as was discussed previously; 
when a = 1 there are additional terms 

£>l , '" l/ "M/*li / <c j73 ( 2 /,0,0,0, / i')^ (/3) (2/ I 0),p (2/))dy. 

The expansion in Lemma 12.61 is valid with \p\ replaced by p and p! replaced by 
1 in view of the following result: 

Lemma 2.7. If f{p) is holomorphic in a sector S which contains the positive real 
axis and if f(p) = \p\ 3 g(p') +o(|/x| 3 |) } then f(p) = g{l)p J +o(|/i| J ). 

Proof. Let C be any closed curve in the sector S. Then for t > 0, 

= t-*- 1 j f(p)dp = t-i- 1 f {\p\ig(p')+o(\p,\i)}dp 
Jtc Jtc 

{\z\ig(z') + o(l)}dz. 

c 

Let t — > oo. By Morera's theorem, |zpf;(z') is holomorphic in 5 so it equals z J g(l), 
the holomorphic extension of its value on {z > 0}. □ 



2.3. The case of a manifold M. The proof of Theorem 12.11 now follows by stan- 
dard arguments [13\ 111) involving a parametrix Pjy on the manifold M constructed 
from the Euclidean parametrices as in Lemma 12.61 In showing that 

(2.q) (P N cf > -(D B +p 2 )- 1 ^p)=0(p- K ) 

for large K, we need to deal with the singularity r~ a in the specific heat. To this 
end, let Rn(v, r, y, s, p) be the kernel of (Db + p 2 )^ 1 — Pn- Then Rm and its first 
derivatives are 0(p~ k ) for large K. Moreover, by construction, the kernel of Pjy is 
zero when r = so the same is true of Rn- It follows that Rm is 0(rp~ K ) for large 
K and Equation pTqf follows. 



3. Heat content asymptotics for Dirichlet boundary conditions 



We adopt the notation of Theorem 11.41 throughout this section. Let B = B-p 
define Dirichlet boundary conditions. We begin by using dimensional analysis to 
express the invariants (3® ^ in terms of a Weyl basis of invariants which is formed 
by contracting indices. 

Lemma 3.1. There exist universal constants e l a so that: 
I dM Pi>% (6 ft A B v)dy = J SM e° a (4, ,p )dy, 

JaM0i™(<l>>P> D > B v) d y = J dM {eKfa, Po ) + el(L aa <p , p ) + el(<p , Pl )} , 
I dM 02%(<t>,P,D.B v )dy = J dM {e i a (cf >2 ,p )+el(L aa ^ 1 ,p )+el(E4> ,po) 
+e 7 a {(j) , pi) + s%(L aa <t>o, pi) +£Q<Ric mm o ,Po) + £a O (L aQ L6b0o, Pa) 
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+E 1 J-{L ab L ab (f> ,p Q ) +4 2 (0 O ;a,PO;a> +£^(^0,^0) + (01 , Pi > }*/ , 

where r is the scalar curvature. 

Proof. Let c > 0. Wc consider a metric g c := c 2 g. We then have for a ^ 1 (see the 
paragraph after (|3.bj) for the case a = 1): 

cte c := c m dx, dy c := c m ~ 1 dy, D c := cr 2 D, 
(3. a) r c := cr, 9 rc = c" 1 ^, 0< )C = c a_i 0i, 

Rcg,c — -^Reg; pi,c — ^ pi • 

We then compute: 

(3((j>,p,D c ,B v )(t)= [ (e- tD ^<f>,p)dx c 
Jm 

= c m f (e c ~ 2tDB &p)dx = c m p(<j),p,D,B v ){c- 2 t). 
Jm 

We take a Z and then continue analytically to the integer values with a^l. The 
interior and boundary terms then decouple so we may conclude: 

OO p 

£ t ( l+ *-«)/2 / £ ^..^^p^DMc^dy 

fc=0 JdM i+j<k 

oo 

fc=0 JajU i+j<fc 

Equating powers oft in the asymptotic series and simplifying yields: 



flfe&jfai, Pi, <r A D, B v )dy = c^-« / Pk™,ij(<f>i>Pj,D, B v )dy . 

Studying relations of this kind is by now quite standard and we refer to [9] for 
further details. We conclude that /3f^f , j is homogeneous of weighted degree A; in the 
jets of 4n and of pj . We now use Weyl's theory of invariants to write down a spanning 
set for the invariants which arise in this way. There is some indeterminacy in these 
invariants as we can always integrate by parts to eliminate tangential divergence 
terms. For example, we have 

(3-b) / (0Q:aa,Po)dy = / (00, pQ:aa)dy = - / (</> :a, P:a)dy . 

JdM JdM JdM 

For this reason, we have eliminated the first two invariants from the formula in 
Lemma 13.11 This completes the proof of Lemma 13.11 for a 7^ 1 . 

If a = 1, the argument is different. In this instance, the regularizing term is 
given by: 

ln(e c ) / ((f>o.c,Po,c}dy c = ln{ce) c(0 o , p )c m ~ 1 dy 

JdM JdM 

= c m {ln(e) + ln(c)} / (fo,Po)dy. 

JdM 

This then yields the modified relation: 

(3.c) lRcg,c(0, p) = c m I Reg (0, p) + ln(c)c m / (<j> , p )dy . 

JdM 
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Recall the notation of Theorem 11.41 One has that: 

00 (—f) n 

/J(0, p, c- 2 D, B v ){t) ~ —T^.M, c- 2n D n p) 



n=0 



fc=0 JdM 

00 - 

+ E E tm / fiZiM^k^'Pi^DiRD^dy 
k=oi+j<k JdM 

00 (—t) n 
' n! 

n=0 

00 . 

+ Vc m - fc t fc / 2 (ln(<)-21n(c)) / p k (<t>, p, D, B v )dy 

k=0 JdM 

OO „ 

fe=0i+i<fc J9M 
Equating terms in the asymptotic expansion then yields: 

(3-d) / flbftjfoi, ft, c~ 2 A Bp)d2/ = c l+J - fe / ftZufapjiD, B v )dy, 

JdM JdM 

t±)l lRcg c ^D n p) = c m \l Rcg (cb,D n p) - 21n(c) / p, £>, Bpjdyl , 



0=-21n(c)/ p 2n+1 (<j),p,D,B v )dy. 

JdM 



The weighted homogeneity of /3f *^ ^ now follows and completes the proof of Lemma 
Theorem O (4) follows from Equations |J5x|) and ([3TH|) □ 



We shall prove Assertions (1) and (2) of Theorem 11.51 by evaluating the normal- 
izing constants in Lemma l3~T1 We begin by establishing some product formulas: 



Lemma 3.2. Suppose that M = Mi x M%, that qm = 9m x +3m 2 > that dM\ = 0, and 
that Dm = Dm ± + Dm 2 where Dmi and Dm 2 are scalar operators of Laplace type 
on M\ and on M2, respectively. Suppose that 4>m = 4>m 1 4>m 2 an d Pm — Pm x Pm 2 
decompose similarly. Then 

(1) (3{<j>M , Pm , D M , B v ){t) = [}{<t>M x , pMi,DMi)(t) ■ f3(4>M 2 , Pai 2 , D M - 2 , BT,)(t). 

(2) J aM P™(<f>M,PM,D M ,Bv)dy 

= E,2n+j=k ^^f- J Ml (fax , (D Ml ) n pM 1 )dx Ml 

x JdM 2 Pj,a 2 (4>M 2 , Pm 2 , Dm 2 , B v )dyM 2 ■ 

(3) The universal constants e l a are dimension free. 

(4) el = 4 )£ i 3 =0 ; ande 1 * = -e° a . 

Proof. Assertion (1) follows from the identity e ~ tDM - B = e ~ tDflf i e~ tDhI ^ B and As- 
sertion (2) follows from Assertion (1). If we take Mi = S 1 , Dm x — —d 2 , (pMi = 1, 
and pM-i = 1, we have that (3(4>m 1 , PMi , Dmi ) (t) = 27r. This then yields the identity 

/ Pk!£(<pM 2 , Pm 2 , D, B v )dy = 2?r / /?f R * 2 {<p M2 , Pm 2 , D M , 2 , B v )dy 2 ■ 

JdM JdM 2 

Assertion (3) now follows. 

We take M 2 = [0, 1] and D 2 = -d 2 . We take 

4>m 2 = Pm 2 = near r = 1, 

Pm 2 = 1 and <$m 2 = r~ a near r — . 
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Since the structures on M 2 are flat, 

f if r = 1 and k > 0, 

Pt M2 (Hh,PM 2 ,D M2 ,B v ){r) = \ if r = and fc > 0, 

[ e° a if r = and fc = 0. 

As the second fundamental form vanishes, the distinction between ';' and ':' disa- 
pears and we may use Eauation (|l.c[) to see that Dip^ 1 = — (p-, aa + Ep). Theorem 
11.41 then implies 

P2{(/>M 1 ,PM 1 ,Dm 1 ) = / {<t>M 1 ,PM 1 -aa+ Ep Ml )d,Xi . 
JdM 

We may therefore use Assertion (2) to see 

p2 A a(4>M,PM-,D M ,B v )dy = e" / (</> Ml , p Jfi;M + Ep Ml )dxi 

3M ' J Mi 

Assertion (4) now follows from this identity. □ 
Next, we evaluate the universal constants e° a . 

Lemma 3.3. 

(1) Ifa^l, then e° a = Tr" 1 /^ 1 "^ (*=s) (o-l)" 1 . 

(2) Lei 7 6e Euler's constant. Then e\ = ^. 

Proof. The proof follows from (12. dp . We note that Assertion (1) also follows for 
1 < a < 2 by the special case calculation in [2]. Assertion (1) then follows for a^l 
by analytic continuation. 

To study the case a — 1 by a special case calculation we let M = [0,oo), let 
= 1 on [0,e] and with compact support in [0, 1), let <j) = r _1 6(r), let p = 1, let 
D = —8%, and let 7 be Euler's constant. As is usual, we work dually and compute 
j3{p, <fi, D, B-p)(t). The halfspace solution of the heat equation with constant initial 
temperature is given by: 



u(r:t) = -7= / e s ds . 
^ Jo 

Consequently, we may compute: 

f3(p,cf ) ,D,B v )(t) = ^ f [ 2Vl e - s \- l Q{r)dsdr 
v Jo Jo 

= \l Hr)d r J 6(r) j ^ e'^ds 1 dr 



if 1 r = 2 r 1 f^Tt 2 

-= / \n(r)Q(r)e--dt -= \ / ln(r)6'(r)dsdr 

Vnt Jo V7T J e J 



= Bi + B 2 , 

where 



B x = --^ij \n(r)e~£dr + J ln(r)(6(r) - l)e"^dr| 
= — — U / ln(r)e 4t + 0(e »* ) 

= __^|y ln(r)e- r2 dr + y m^K^r j + 0(e - *r) 



= -§ ln(i) + j- ln(2) - ^ y hi(r)e- r dr j + 0(e~ w) 
ln(r)9'(r) j y e^'ds - J ^ e'^dsj dr 
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]n(r)9'(r)dr + 0(e~&) 



= -]u(r)0(r) +J r -1 8(r)dr + 0(e _ *f) 

= ln(e) + y r- 1 6(r)rfr + 0(e-fe). 

This then yields the expression 
(3(p,<t>,D,Bv)(t) 

/■CO pi 2 

= iln(^)-ln(2)- ^y ln(s) e - s2 ds + J r _1 e(r)dr + 0(e - w). 

Since <^> is compactly supported in [0, 1), the heat content for corresponding problem 
on the interval [0,1] is the same as for [0, oo) up to 0(t°°). Assertion (2) now 
follows. □ 

We continue our study by index shifting: 
Lemma 3.4. 

(1) Assume Vg r <& = on C € . Set (j) := x(r)$(y)r io_7 . Then 

Iom Pk%,<ojfr ft"' A B v) d y = S dM /%,7-*>,o,i(*> Pi' A A>)dy • 

(2) I/a ^ 1, iften 4 = 4_ x , 4 = e° a _ 2 , e 5 a = and 4 4 = 4_ t . 

(3) ej = ef = -I, e\ = \, and e} 4 = 0. 

Proof. Set := x( r ) < E > (y) r4 ° By Theorem II .41 with a — 7 and with a = 7 — io, 



>+j<k JdM 

JdM 



We set k = £ + io and equate powers of i to establish Assertion (1); Assertion (2) 
then follows from Lemma 13.11 and Assertion (3) from Theorem 11.31 □ 

We continue our study with a functorial property that exploits the fact that we 
are working in a very general context; we are no longer working with the scalar 
Laplacian! Even if one were only interested in the scalar Laplacian, it would be 
necessary to consider general operators of Laplace type in order to use this functorial 
property! Let O = 1 on [0, h] and with compact support in [0, 1). 

Lemma 3.5. Let T" 1_1 be the torus with periodic parameters (yi, y m -i). Let 
M = T m ~ l x [0, 1]. Let f a G C°°([0, 1]) satisfy f a (0) = and f a = near r = 1. 
Let S a G K. Set 

ds 2 M =Y,a e2fa(r)d Va°dy a + drodr 1 p := e~Z a /.M, 
Dm ■= - E a e- 2fa{r) (d 2 ya + 5 a d y J - d 2 r , <i> := 6(r)r- Q . 

(1) If k > 0, then J dM /3g£(0, p, D M ,B v )dy = 0. 

(2) -|4-4 = o. 

(3) -K4 + 4 2 ) = o. 

C4> _I e 4 + I £ 6 _ i £ 7 _ 9 = q 

V^V 4 a T 2 a 4 a c a u - 

(- ip4 , 1^5 , 1_6 , 1-7 , 1_8 i _10 _ n 

V^V 8 c a^2 a^4 G a'8 a 2 a^ c a 

(6) -4 + 41=0. 
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Proof. We use — <9 2 on [0,1] and Dm on M. Since 4> vanishes near r = 1, this 
boundary component plays no role. Let !tg(r; i) be the solution of the heat equation 
on [0, 1] with Dirichlet boundary conditions. Since the problem decouples, UB(r;t) 
is also the solution of the heat equation on M with Dirichlet boundary conditions. 
The Riemannian measure 

dx = \J det gijdydr — e^ a ^ a dydr . 

As p — e~^a fa , pdx = dydr. We suppose a^l. Since vol(T m_1 ) = (27r)" l_1 , 

l3{0,p,D,B v )(t)= J u B (r;t)pdx ^ {2n) m - 1 J u B (r;t)dr 

= (2^) m " 1 /3(0, 1, -3 2 r ){t) = (2 7 r)'"- 1 (X Rcg (0, 1) + e° a ) + 0(t n ) 

for any n since the structures are flat on the interval. Note that Dp = 0. Assertion 
(1) now follows. If a = 1, the computation must be modified to take the In term 
into account; this does not affect the computation of /3f M for k > 1 and the desired 
conclusion follows similarly 

To apply Assertion (1), we must determine the relevant tensors. The formalism 
of Lemma ll. II is crucial at this point as the connection defined by the operator Dm 
is no longer flat. We have: 

= 5<5a, &a = ~^>a = — 5^°' 

w «> = ~\ Ea /a> &m = = \ Ea fa ■ 

We compute: 

fa = {v 9r (r^)}U = {{d r - |E /:)(i)}|9M - 4Ej:(o), 
(3e) = |(E Q /^(o)) 2 -iE a /:(o), 

Po = 1, 

Pi = {V 9r (p)}U = {(dr + § Ea f'a)(e- E » fa )}\dM = ~\ Ea /a(0), 
P2 = ±{(^)MU = \{(dr + \ Ea fafi^ E " U )}\bM 

= |(E a /^(o)) 2 -iE Q /:(o). 

It is straightforward to compute that we have the following relations when r = 0; 
we refer to Lemma 2.3.7 9 for further details: 



(3.f) 



E 2^oA ^ 4C0/J |£oi' ^aa Ea ■) "a J 

RiC mm = - Ea((/a) 2 + /a ) > L aa L bb = (E a f' a ) 2 i L ab L ab = Ea(fa) 2 



Considering the term Ea fa m /^laf yields Assertion (2), considering the term 
E Q <5 2 in /3f cf yields Assertion (3), considering the term E a fa m /^fcf yields Asser- 
tion (4), considering the term (E a f a ) 2 ln 02 of yields Assertion (5), and considering 
the term Ea(fa) 2 in 02% yields Assertion (6). □ 

We continue our discussion: 

Lemma 3.6. 

(1) 1/^0 = 0, thend t f3(cj ) ,p,D,B v )(t) = ~(3{^Dp,D,B v ){t). 

(2) 4 = 0. 

(3) //a ^ 1, 4 = 3^4, 4 = -3^4, and 4 4 = 0. 

(4) = 2e° + 1, ef = -e? - |, and e^ 4 = 0. 

Proof. Assume pa — 0. By Equation Ql.cp we have 

(3.g) - (L»p)o = 2/? 2 - iaaPl ■ 
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We compute that: 

d t p{4>,p,D,Bv){t) = -{De~ tDB cl>,p) = -{e- tDB 4>,Dp) 
= -P(<j>,Dp,D,B v )(t), 

where the middle equality is justified as po = 0. This proves Assertion (1). We use 
Assertion (1) to see 



f (3™(^p,D,B v )dy = Q. 

JdM 



IBM 

Since there is no restriction on px, we conclude that e 3 = which establishes 
Assertion (2). Furthermore, if a ^ 1, we may conclude 

( 3 . h ) i±|=a f pMfa /0) Dj Bv )dy = - [ (3 d k \ a {4>, Dp, D, B v )dy . 

JdM JdM 

We set k = 2 to see that 

e° Q (fo,(Dp)o) 

dM 

= T 2 / {{£l(<Po,P2)+e s Q (L aa fo,p 1 )+e 1 a 4 {fa,p 1 )}dy. 

JdM 

Assertion (3) now follows from Equation ( |3.g[ ). If a — 1, we have: 

\±{(f> ,(Dp) } +el(<po,p 2 ) +el(L aa (j) Q ,p 1 } + e 1 4 (0i, pi)} dy . 

IdM v > 

Assertion (4) now follows. □ 

Assertions (1) and (2) of Theorem 11.51 will follow from Lemma I3T21 and from the 
following result: 

Lemma 3.7. 

(1) Suppose that a =/= 1. Then: 

(a) e a = c a . 

(b) e\ = Ca-x, e 2 a = -\c a -i, and e 3 a = 0. 

(c) e% = c Q _ 2 and e b a = - ±c Q _ 2 - 

IA\ c 6 a— 3 c 7 2 j P 8 1 „ 

\ u ) fc a — 2(a-l)(a-2) L "- 2 ' t a ~ ( Q _ 1) (a-2) Ca ~ 2 uaut a — (a-l)(a-2) Ca ~ 2 ' 

( e ) £ a - ~ 4(a-2) C "- 2 £ " _ 8(a-2) C °- 2 - 

(f) £l a = ~ 4(a-2) C "- 2 ' ^ = 2(a-"7(a-2) Ca ~ 2 ' e af = ' flrlrf ^ = ' 

(2) Let a = 1. T/ien: 

(a) £ ? = | 7 . 

(b) e\ = —fa, e\ = ^, and el = 0. 

(c) 4 = -l, s\ = I. 

(d) e\ = |7, e\ = 7 + 1, and e? = —^7 - |. 

(e) e? = and ej° = 0. 

(f) e} 1 = 0, ef = ±7, e} 3 = 0, and e} 4 = 0. 

Proof. Let a 7^ 1. Assertions (la) and (2a) follow from Lemma l3~3l Assertions (lb) 
and (2b) follow from Lemmas l3. 41 13.51 and 13. 61 Assertions (lc) and (2c) follow from 
Assertion (lb) and from Lemma [3~4l Because sT(s) — T(s + 1), we have: 

a — 3 

rCa-2 ■ 



2(a-l)(a-2) 
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Assertions (Id) and (2d) now follow from Lemmas 13.21 and 13.61 We use Lemma 
to establish Assertions (le) and (2e) by computing for a / 1 that 



£ 9 = -±£ 4 + ±£ 6 - ±£ 7 



„ f_l _ 1 a-3 1 2 \ _ q-1 

t-o-2V 4 4 (a-l)(a-2) 4 (a-l)(a-2) J ~~ 4(q-2) Cq - 2 ' 



=- 10 = _Ip4 _ 1„5 _ 1.-6 _ 1_7 _ 1, 
fc a 8 Q 2 a 4 a 8 a 2 



_ _ / 1 i 1 i 1 a-3 1 1 i 1 1 \ _ a-1 - 

— <--a-2\ 8 -I- 4 T- g ( Q _i)( Q _2) 4 (a-l)(a-2) 2 (a-l)(a-2)/ — 8(a-2) L "- 2 

and for a = 1 that 

= -it" 1 ) + 5(h) -i(7 + l)=0, 

4° = -i(-i) - - k(h) - i(7 + 1) - - 1) = o- 

Assertions (If) and (2f) follow from Lemmas 13.51 13. 2[ and 13.61 □ 

4. Heat content asymptotics for Robin boundary conditions 

Section [4] is devoted to the proof of Assertion (3) of Theorem 11.51 Let B — B-ji 
define Robin boundary conditions. We clear the previous notation concerning the 
constants e l a . Recall B-np = pi + Spo. Lemma 13.11 extends immediately to this 
setting, after including the additional tensor S in the Weyl calculus, to yield: 

Lemma 4.1. There exist universal constants e l a and <P a so that: 

(1) I 9M P, D, B n )dy = J 9M e° a (fo, Po )dy. 

( 2 ) /aAfC^ft- -^)* = IoaA^uPo) + eKLaafo, po) + 4(</> ,Pi) 

+d 1 a {<p ,B lzP )}dy. 

(3) f dM P$M(<j>,p,D,B n )dy = f aM {e i a {fc,Po)+el{Laa<t>i,Po)+e%{E<p ,p ) 

+E 7 a (<fi , p 2 ) + el(L aa (f) , pi) + £ 9 l (Ric mm o ,Po) + e]^{L aa L bb (j) , p ) 

+£ 1 a 1 (LabLab<f>0,Po) + £ a 2 {4>0;a-, P0-a) + S^ 3 (rfo , p ) + ^^(^l: Pi) 

+ (d 2 a cj) 1 + dlScfx) + d 4 a L aa (j> a , B nP )}dy. 
We begin our analysis by showing that all the constants e l a vanish: 
Lemma 4.2. 

(1) IfB nP = Q, thend t f3{^,p,D,Bn){t) = ~P{^Dp,D,B n ){t). 

(2) 4 = 4 = 4 = 0. 

(o\ P 4 _ p 5 _ p 6 _ P 7 _ p 9 -10 _ -11 _ -12 -13 _ n 
ct t a fc a fc o fc a t a fc « fc a U- 

(4) 4 = 4 = 44 = 0. 



Proof. Assertion (1) follows using the same arguments used to prove Lemma 13.6 
(1); Equation (|3.hp then generalizes to become 



/ $%(<t>,P,D,B n )dy = - [ p a k ™ 2a (^Dp,D,B n )dy. 

JdM JdM 



We take S = and p\ = 0; po and p2 are then arbitrary. Since fi^f a — an d 
/3£M q = 0j Assertion (2) follows. This implies that /3$™ (<f>, p, D, B n )' = and a 
similar argument now establishes Assertion (3). We now take S = — 1 and po — 
pi = 1 to establish Assertion (4). □ 

The constants d^, d\, and d\ can be determined by a 1-dimensional calculation. 
We adopt the following notational conventions. Let M := [0, 1], let A := d x + b 
where 6 € C°°(M) is real valued, let A* := -d x + b, let D x := A* A, let D 2 := AA* , 
and let Bn4> '■= A4>\qm- The inward unit normal is d x near x = and — near 
x = 1. Thus this is a Robin boundary condition with 5(0) = 6(0) and 5(1) = —6(1). 

Lemma 4.3. Let a G C — Z iwt/i Re(a) < 0. Adopt the notation established above. 
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(1) d t (3(<j>,p,D x ,B n )(t) = -P(A^,Ap,D 2 ,Bx,)(t). 

(2) J BM P d u%{4>, P, D u B n )dy = J dM /3£^ Q+1 (A0, Ap, D 2 ,B v )dy. 

(3) di- ~ • 



2-c* 

2 



(4) dlfa + 450o = - 3^c Q {(l - a)0! + S<f>o}- 

Proof. We generalize the proof of Lemma 2.1.15 [9] where a similar result is es- 
tablished for a = 0. One has that Ae~ tDl - B ^ = e~ 2 ' B v A on sufficiently smooth 
functions. Thus we may establish Assertion (1) by noting: 

d t {e- tD ^4>,p) = -{A*Ae- tDl < B -K-(j),p) = -{A*e- tD2 > B vA4>,p) 



= - e 



2 - B v A(j>,Ap) 



where the middle equality is justified by the boundary condition B-p. 
Since r a cf> G C°°(V), we have r a+1 A<f> G C°°(V). Thus 

oo - 

~ -£^ (1+£ - (Q+1))/2 / fl% +1 {M,MD,ih.)du. 

Setting ^ = fc — 1 and equating coefficients of ^-i-fe-a)/ 2 yields Assertion (2). 
The operators 

^ 1 = -(^ + (&'-6 2 )) and D 2 = -(^ + (-{/- 6 s )) 

determine flat connections. We suppose that 6, 0, and p vanish identically near 
r = 1 so only the point r = is relevant. We set 0_i := and expand: 

oo oo 

^ r ^ a ' A ^ ~ £ {(* - a ^ + ri ~ a ~ l ■ 

i=0 i=0 

It now follows that (A0)o = — a0o an d (/4-0)l — (1 — a)0i +&0o- We apply Assertion 
(2) with k = 1 and k — 2 and we apply Theorem ll.5l (1) to see: 

{d 1 a <t>o,Ap)dy = c a+x / (-a^> ,Ap)dy, 

dM £ — OL Jqm 

(401 + c^0o, = - » 2 c a / ((1 - a) 0i + 60o, Ap)dy . 

dM 6 — a J dM 

Assertions (2) and (3) now follow. □ 



We extend Lemma 13.51 to the setting at hand to complete the proof of Theorem 
0(3). 

Lemma 4.4. Adopt the notation established in Lemma VJ . 5\ Let S = -z ^2 a f' a define 
Robin boundary conditions. Then 

(1) f eM Pm-Dm, Bn)dy = 0. 

(2) di = -3^c a . 

Proof. Taking into account the change in the connection, we have that B-jz on M 
agrees with the pure Neumann operator Bj^f defined by S — on [0, 1]. Since one 
has that pdx = rx{r)dydr, 

/3(0, p, D M ,B K )(t) = (2 7 r) m - 1 /3(0, rx, -d 2 rl B N ){t) . 

Assertion (1) now follows as f g r Q ^ /?f.* / (0i r X> —d^,Bj^f)dr = 0. 

To prove Assertion (2), it is simply a matter of disentangling everything. We use 
the equations of structure derived in the proof of Lemma 13.41 to see: 

0o = l, ^i = -iE a /a! Po = 0, pi = 1, 

S = 2 So /a, ^«» = — /Ca fa- 
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We may now compute: 

0= / {-&C a {(l- a )(-±Zaf'a) + hZati)+<t i a (-Zaf'a)}dV- 
JdM 

It now follows that d A a = — -^z^c a . □ 
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